%求解控制序列为a_min时候的相关参数
%t_f = x(1)
F = @(x)[12 - 4 * x(1) - 12 * (1 - 0.15 * sin(0.2 * x(1)));
    30 + 12 * x(1) - 0.5 * 4 * x(1) * x(1) - (2 * 10 - 2 * 1.8 * 12 * (1 - 0.15 * sin(0.2 * x(1))) + (100 + 12 * x(1) + 10 * cos(0.2 * x(1)) - 10))];

%创建初始点 x0 
x0 = [0];
%设置选项以返回迭代显示。
options = optimoptions('fsolve','Display','iter');